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Information on the decay process of nuclides in the superheavy region is critical in investigating new elements 
beyond oganesson and the island of stability. This paper presents the application of a random forest algorithm to 
examine the competition among different decay modes in the superheavy region, including a decay, 87 decay, 
BT decay, electron capture and spontaneous fission. The observed half-lives and dominant decay mode are well 
reproduced. The dominant decay mode of 96.9% of the nuclei beyond ?'?Po is correctly obtained. Further, 
a decay is predicted to be the dominant decay mode for isotopes in new elements Z = 119—122, except for 
spontaneous fission in certain even-even elements owing to the increased Coulomb repulsion and odd-even 
effect. The predicted half-lives demonstrate the existence of a long-lived spontaneous fission island southwest 
of ?95F] caused by the competition between the fission barrier and Coulomb repulsion. A better understanding 
of spontaneous fission, particularly beyond ?99FI, is crucial in the search for new elements and the island of 
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stability. 
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I. INTRODUCTION 


Limitations of the nuclear landscape [1, 2] have always 
s been an intriguing topic. Exotic nuclear properties, for ex- 
ample, the shell evolution [3—6], 4n resonant state [7, 8], and 
4p unbound state [9], emerge at the limits of nuclear stability. 
The discovery of new elements (nuclides) involves the follow- 
ing three problems: production, separation, and identification 
[10]. Because the nuclei are unstable and have relatively short 
half-lives, appropriate probes must be utilized. Characteristic 
decay modes [10, 11] are commonly employed as a probe to 
signal the existence of exotic nuclei. Therefore, investigating 
and predicting the dominant decay modes of the unknown nu- 
clides is crucial. The nuclear binding energy and half-life are 
key parameters for determining the decay mode of a nucleus. 
The former measures the stability of nuclides by using energy 
criteria, and the latter describes the possibility of decay. 

Both microscopic and macroscopic methods have been 
used to study the nuclear binding energy [12-16] and partial 
half-life of each decay channel, including a decay [17-19], 8 
decay [20, 21], spontaneous fission [22, 23], protons emission 
[25] and neutron emission [24], etc. Microscopic theories be- 
gin with nucleon-nucleon interactions, which can be based 
on either realistic or phenomenological models. The macro- 
scopic theory uses selected variables with physical consider- 
ations to construct semi-empirical formulas and fit the exper- 
imental data, and it entails the risk of overfitting and inappro- 
priate parameters. In addition, exotic nuclei may significantly 
deviate from the general fitting and be identified as outliers. 
Decreasing the deviation between theoretical predictions and 
the observed results remains a critical issue. 

With advances in computing and storage, efficacious ma- 
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chine learning algorithms with diverse applications have been 
proposed [26, 27], e.g., nuclear properties [28-30], fission 
yields [31-35], spectra decomposition [36], radiation effect 
[37], neutrino experiment [38], and other nuclear techniques 
[39-42]. As summarized in a recent colloquium, estimat- 
ing the residuals of nuclear properties using machine learn- 
ing algorithms is a powerful strategy [43]. A neural net- 
work has been used to compensate for the residuals of nu- 
clear masses [44—47] and nuclear charge radii [48—50]; this 
has been achieved through structural optimization and care- 
ful selection of the input parameters with definite physical 
interpretations. The applicability of the decision tree (DT) 
has been verified via training and testing with residuals of 
the binding energies [51]. However, random forest (RF) [52] 
algorithm, developed from the DT algorithm, has not been 
tested for determining the nuclear mass or the partial half- 
life of a specific decay channel; in this regard, semi-empirical 
formulas have suggested several major components but with 
residuals. Machine learning algorithms can include possible 
features to realize the training for residuals, whereas RF, with 
bootstrap sampling, not only avoids overfitting but also con- 
siders the correlation between the data combinations and sev- 
eral features. Thus, RF exhibits increased robustness and is 
conducive to extrapolation. The amount of computation in- 
creasing in accordance with the number of trees in the forest 
and size of the dataset, as well as the difficulty of model in- 
terpretability, may limit its applications. 

This study entailed the application of an RF machine learn- 
ing algorithm to analyze the major decay modes of heavy and 
superheavy nuclei. The competition between a decay, 8 de- 
cay, and spontaneous fission (SF) of new elements, as well as 
the possible long-lived island in the superheavy region were 
examined. 


Il. METHOD 


This study focused on the Z > 84 and N > 128 regions. 
The partial half-lives of the o decay, 8~ decay, B+ decay, 
electron capture (EC), and SF were calculated using semi- 
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empirical formulas, and the residuals of each formula were 
then trained using the RF algorithm. The minimum partial 
half-life of a mode corresponds to the dominant decay mode. 


A. Decay Half-life Formulas 
The universal decay law (UDL) [53, 54], 


log igT1/2,0 —aZ4(Z = Za) V u/Qa 
Y b uZa(Z m ZA? RE (A m Api) 


+c, 


() 
is used to fit the a-decay half-life. Za, Aa, Qa, and u = 
A,(A — A,,)/A denote the proton number, mass number of 
the o particle, a decay energy, and reduced mass, respec- 
tively. The channel is assumed to move from the ground state 
to the ground state. 

A three-parameter formula (denoted as SF3) was used for 
the SF as follows: 


logio Tsp = à a 


which was proposed based on several existing formulas 
[22, 23, 55-58], where v represents the blocking effect from 
unpaired nucleons; its value is 0 for even-even nuclei and 2 
for other nuclei [55]. « has a value of 2.6 [22, 59], I = 444, 
and a, b, and c are the fitting coefficients. Eq (2) is separately 
fitted to nuclei with Z « 104 and the remaining because of a 
systematic difference, as shown in TABLE 1. T! sp of the 
nuclei with Z « 104 increases significantly with a decrease 
in Z because the Coulomb repulsion decreases. The relatively 
long Ts (> 10? s) of certain nuclei in this region cannot be 
universally described currently and were not used to fit Eq.(2) 
because the competition for the SF is significantly weak com- 
pared to other decay modes. 

The 6 decay half-life was estimated using the formula 
given in Refs. [21, 60]. Assuming that the ground state 6 
decay is an effective Gamow- Teller (GT) transition, the par- 
tial half-life can be expressed as follows: 


logioT1/25,5 = logig&1 — logigfo — log, ) Ber, (3) 


EN 22? h* 1n 2 E 


where &1 më? = 6147 s, fo is the phase-space factor, 


and Bor is the GT-reduced transition probability [60]. As 
regards EC, the phase-space factor is deduced as follows: 


EC 3 2 2 4 
whereas for the G~ decay, it is 
BE F(E — 10E + 15E — 6)2«(Z F 1)/137 (5) 
o m 


30(1 — exp(E22(Z 7 1)/137)) ? 


where Eo is the renormalized 5-decay energy. Because Qg 
provided by AME2020 [61] is the difference in the atomic 
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mass, the electron mass should be reconsidered as follows: 


Q g^ + Mec? 
E, = 
0.8 Mec? 
Qg- T mec 
Eag- = 6 
ba a3 (6) 
2 
— mec 
Pie Qzc i 
mec 
Finally, logigBcr is estimated às the average 
logio(foT1/2,8/K1). The fitting results are listed in 
TABLE 1. 


B. Random Forest Method 


RF is a fusion of the DT and bootstrap algorithms. DT 
is a non-parametric supervised learning algorithm. For a 
dataset consisting of S samples of J features (variables) 
((01,...,07)5, s € [1,S]} and object (observable) (ys, s € 
(1, S]}, it establishes a binary tree structure that divides the 
dataset into L subsets based on the values of the features; 
each subset is called a leaf. This partition seeks to minimize 
the RMSE 


S 
RMSE = So (ys — f (O15 97)? (7) 
g=] 


Ul = 


of the entire dataset by assigning a value to each leaf. 

Bootstrap is a statistical method based on the concept of 
random resampling with replacement, through which possi- 
ble combinations and weights of data are automatically con- 
sidered [62, 63]. Each time a new dataset is obtained, a new 
DT is trained and used to predict the object of each sample in 
the entire dataset. By repeating this process M times, a forest 
of M trees is obtained. The final predicted value of the object 
for a sample is the average of the results calculated by all the 
trees in the forest. Because each tree is trained by part of the 
samples in the dataset, the value for each sample predicted by 
the forest is an average of the interpolation and extrapolation; 
this decreases the divergence when the calculation is imple- 
mented for the unmeasured nuclei. The open-source Python 
library scikit-learn [64] was used for machine learning. The 
forest was assumed to be composed of 10° trees so as to de- 
crease the dispersion of the RMSE in this study. 


IIl. RESULTS AND DISCUSSION 


The residuals of the decay formulas of a decay, 67 decay, 
Bt decay, and EC were trained using the RF with features Z, 
N, A, the odevity of Z and N, and the decay energy. Be- 
cause the decay energy cannot be defined for the SF, the fis- 
sion barrier (FB) obtained from Ref. [65] was used to replace 
the decay energy in the feature set to consider the deforma- 
tion effect. The number of leaves chosen for this study was 
11, which was the same as that used for training the bind- 
ing energy in this region in our previous study [66]. Fig. 1 
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Table 1. Coefficients and corresponding root-mean-square error (RMSE) of the UDL, SF3, and Eq. (3) when fitted to the nuclei with Z > 84 
& N 2 128. The RMSEs of the RF-trained UDL, SF3, and Eq. (3) are listed in the last two columns. The Weizsácker—Skyrme (WS4) [12] 
and universal nuclear energy density functional (UNEDFO) [13] in the subscript indicate the sources of predicted energies. 


a (log1o Ber) b c RMSE RMSEnr,ws4 RMSEnr,uxEDrFo 
UDL 0.407 -0.382 -23.896 0.883 0.598 0.669 
SF3 (Z « 104) -1.129 -6997.113 79.803 3.070 1.195 1.195 
SF3 (Z > 104) -1.363 -13272.729 113.415 1.267 0.825 0.825 
Eq. (3)e+ 1.378 - - 1.957 0.439 0.437 
Eq. (3)g- -1.819 - - 1.451 0.656 0.667 
Eq. (3)gc -2.112 - - 2.360 0.971 0.996 


compares the residuals of these decay formulas before and af- 
ter RF training. Two conditions were assumed to determine 
the outliers: 1) located outside the dashed line with the cor- 
responding color, which indicates that the scatter deviates by 
twice the RMSE from the experimental log; 5; and 2) the 
[logio (T1/2,ca1/T1/2,exp)| value is larger than 3, which indi- 
cates that the calculated value is three times that of the mag- 
nitude of the experimental value. Thus, missing (adding) the 
outliers owing to the significantly large (small) RMSE can be 
avoided. After training, the biases of the outliers of these de- 
cay formulas were significantly reduced, and the RMSE of 
the formulas decreased (TABLE 1), as expected. The condi- 
tion of the outlier is not too strict because the aim was not 
to maximally decrease the RMSE but to reach an appropri- 
ate scale wherein the dominant decay mode can be described. 
The same features and number of leaves of the RF were cho- 
sen in this study to train the residuals of the different decay 
formulas; this avoids overfitting while seeking an extremely 
small RMSE. 

In total, 445 nuclides with measured partial half-lives and 
branch ratios of the five decay modes were obtained from 
NUBASE2020 [67]. The dominant decay modes and partial 
half-lives of the nuclides are illustrated in Fig. 2(a,b). A long- 
lived a decay valley from 32°Raj38 to 22. Cf153 lies between 
the narrow 6+/EC decay band and the neutron-rich 67 re- 
gion. The half-life of the nucleus decreases with increasing 
distance from this valley. The southwest was dominated by 
a decay, whereas the southeast was dominated by 9^ decay. 
In the northwest, 8* decay and EC compete with a decay 
and lose after Z increases. In the northeast, a decay and SF 
compete with one other, and the region extending from the a 
valley appears to be dominated by the SF. Although the dis- 
183 tribution of the dominant decay modes demonstrates a clear 
184 boundary, the minimum partial half-life was smooth. 

Among the 445 nuclides considered, 341 (104) nuclides 
186 had known (unknown) corresponding decay energies. The 
15; nuclides with unmeasured masses were calculated using WS4 
188 [12] and UNEDFO [13] to estimate the partial half-lives. The 
189 results of the RF are presented in Fig. 2(c-f). The calcu- 
190 lated results sufficiently agree with the experimental results 
191 because the dominant decay mode is correctly described for 
192 431 and 427 (96.9% and 96.0%) nuclei when the RMSE of 
193 log19T1/2 of the dominant decay mode is 0.62 and 0.67, re- 
194 spectively. Nuclides, for which the dominant decay mode was 
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Fig. 1. (Color online) Comparison of the residuals of UDL (a de- 
cay), SF3 (SF), Eq. (3) (8^, 8* and EC) before and after RF train- 
ing. The dashed lines denote twice the RMSE of the corresponding 
formulas. 


inconsistently described, generally have two competitive de- 
cay modes. For example, the œ and SF branch ratios of ?°°Rf, 
?6? Dp. and 78°F] were approximately 50%. Meanwhile, the 
liquid drop model trained by RF [66] was also applied to ob- 
tain the energies; the model afforded consistent results that 
are not presented herein. 


The N — Z — 48 chain, where the partial half-lives of the 
a decay, 8 decay, and SF are comparable, was used for a spe- 
cific comparison. The results of the density-dependent clus- 
ter model within the anisotropic deformation-dependent sur- 
face diffuseness [68], Royer formula [70], modified Swiate- 
2% Cki’s formula [22], nuclear liquid drop model [23], and de- 
207 formed self-consistent Hartree—Fock mean-field with Skyrme 
28 forces and pairing correlations [70] were compared with the 
209 experimental values; the results are presented in TABLE 2. 
2:0 The models afforded similar results and deduced a consistent 
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(Color online) Dominant decay mode (left panels) and minimum partial half-lives (right panels) of the œ decay, 87 decay, 8+ 


decay, EC, and SF. (a,b) Experimental data in NUBASE2020. (c-f) The results predicted through RF; WS4 [12] and UNEDFO [13] denote 
the sources of predicted energies. Specifically, the FB is used to replace the decay energy to learn SF. The nuclides, for which the predicted 


partial half-life is longer than 10^ s, are marked by a star. 


dominant decay mode. 

The accuracy of the obtained energy is crucial for half-life 
calculations. If the calculated energies of WS4 and UNEDFO 
replace all the experimental energies, the number with a con- 
sistent decay mode compared to the experiment reduces to 
72.6% and 64%, respectively, and the RMSE of logioT! /» 
increases to 2.07 and 2.64. The difference in the results ob- 
tained using the energies of the two models is owing to the 
accuracy because the RMSE of the mass of WS4 is approx- 
imately 0.3 MeV [12], whereas that of UNEDFO is approxi- 
mately 1.45 MeV [13]. This also leads to differences during 
extrapolation. The consistency rate of the dominant decay 
mode between the energies calculated using these two mod- 
els decreases from 82.2% to 66.2%. More accurate and pre- 
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cise measurements of the decay energy will aid in theoretical 


2: predictions. In addition, WS4 and UNEDFO may lose their 
227 predictive power after training with machine learning. Train- 
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ing the WS4 and UNEDFO binding energies with features 
Z, N, 6, and P, which sufficiently describe the residuals in 
Ref. [45], improves the energy description but decreases the 
consistency in the dominant decay mode by several percent, 
which is considerable compared to the 23.4% rate of the the- 
oretical energies among all nuclides (104/445). 

The SF is important for investigating the half-lives of su- 
perheavy nuclei. As shown in Fig. 2(c, e), the dominant 
decay mode of the unknown nuclides is determined in ac- 
cordance with the competition between the SF, o decay, and 
B~ decay. The major competition is between the SF and 67 
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Table 2. Comparison of the experimental partial half-lives of the N — Z = 48 chain with the values calculated by different models. 


[23] 


(22) (70) 


Nucl IgTa exp lgTa cai gi 1T IgTsr exp lgTsr cat ZTE?! IgTE?) IgTo4 exp TEGexp laeva TR ET 
7TUMCf 3.193 3298 3.000 3334 - - - - - 3.670 - 4.426 3.403 
246Es 3.658 4.464 - - - - - - 2.698 - 2.213 - - 

245 Em 1.538 — 1.637 3.358 1.73] 4.538 4.678 4.739 5.069 - - - - 2.025 
250Md 2.887 2.567 - - - - - - 1.764 - 1.809 - - 
?5?No 0.562 0.541 0.253 0.675 0.897 1.799 1.499 2.119 2.351 - 2.223 - 1.822 
2541r 1224 1.159 - - - - - - 1.627 - 1.577 - - 
256RF 0.328 0.082 -0.198 0.250 -2.179 -1.382 -1.071 0.519 - - - - 1.898 
258Db 0.530 0.037 - - - - - - 0.780 - 1.349 - - 
260Sg -1.768 -1.940 -2.300 - -2.157 -2.862 -2.811 -2.301 $ = 2 = - 
decay for neutron-rich nuclides, and between the SF and a 

decay for neutron-deficient nuclides. Existing experimental FB(No T1/2,sr,Exp) 

data demonstrates a long-lived a decay region from 229 Ra13s FB(BR 2 1%) 

251 n u s FB(BR < 1%) 

to 3 Cf155, lying between the 6* and £^ decay regions, logioTı/2 sr eap (BR > 1%) 

and ending with the SF. The proposed models correctly de- log1oT1/2,5F,exp(BR < 1%) 

scribe this phenomenon. In the long-lived region, after N 

exceeds 154, the blue band shown in Fig. 2(d, f) indicates " 
half-lives of approximately 107—107 s. At the southwest cor- z e 
ner of Z = 114 and N = 184, nuclides in the circle have Z na 
a longer half-life than those in the surrounding area. This is cà E 
because of the high FB in this region, which leads to a longer E E 


Tijya,sr. Fig. 3 compares the evolution of FB and measured 
T\/2,sr along the mass number. The FB decreases with A 
before A — 230 and subsequently behaves as a sinusoidal 
wave oscillating between 2 and 10 MeV. Apparently, an FB 
threshold exists, below which SF occurs. Nuclides with rel- 
atively long T1/5 sr. generally have small SF branch ratios. 
In addition, the FB of nuclides with SF branch ratios less 
than 1% were mostly higher than those with SF branch ratios 
greater than 1%, which implies that the higher the FB, the 
weaker the SF. However, considering only the nuclides with 
SF branch ratios of less than 1% or greater than 1%, the cor- 
respondence between FB and T/2 sr becomes significantly 
more complex. 


The nuclides with partial half-lives predicted to be longer 
than 104 s are marked with stars in Figs. 2(d, f), which sug- 
gests 250.252,254 yy, 260,261 p. 261-264 14. and 265] y for fu- 
ture measurements. No experimental value of the half-life of 
250Cm was suggested in NUBASE2020 and was thus extrap- 
olated in this study. In the NNDC, SF was shown to be the 
dominant decay mode, and its half-life was recommended to 
be 8300 years, which is relatively long. Although the calcu- 
lations in this study underestimate the NNDC value, the long 
half-life and dominant decay mode are reproduced. In addi- 
tion, the upper limit of the half-life of ??Cm was proposed to 
be two days in Ref. [69], which was not updated since then 
(1966), whereas the current study estimates a half-life of 1.43 
days. No experimental half-lives were previously reported for 
260,261 ps. 261-264Md, and 26°Lr. However, their nearby iso- 
topes have long half-lives, such as ??"Es (7.7 days), ?9?Md 
(31.8 days), 2°°Md (1.6 h), 2°°Md (51.5 days), ??" Md (5.52 
h), and ?96] (11 h). Moreover, the Es, Md, and Lr isotopes 


28 


282 
283 
284 
285 
286 
287 
288 
289 
290 
29 


292 
293 
294 
295 


2 


© 


6 


iir T. p cic 


210 220 230 240 250 260 270 280 290 300 
A 


Fig. 3. (Color online) Evolution of T; /2,sp and FB along the mass 
number. The datasets are divided according to whether the corre- 
sponding 7:/5»,sr is measured and whether the branch ratio (BR) of 
SF is less than 1%. 


are located in the extension of the narrow long-lived region 
from ??6Ra to ?9! Cf, substantiating that the Es, Md, and Lr 
isotopes are candidates with long partial half-lives. Obtain- 
ing more measurements is also suggested; for example, the 
?52Cm data has not been updated for more than 50 years. 

A comparison of all the possible decay channels is lim- 
ited by the accurate description of each channel and the ob- 
served data; note, the SF mechanism remains unclear, such 
as its dependence on the FB or deformation. The effect of 
the quadrupole deformation parameter (€2) [15] on the half- 
life estimation was then investigated. If the FB is replaced 
with £2 during RF training, more nuclides are predicted to 
have longer half-lives. Further investigations should be con- 
ducted to understand the dominant factors that contribute to 
the half-life of SE. The FB combines the contributions of mul- 
tipole deformations and thus presents a stronger quantum ef- 
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Fig. 4. (Color online) Comparison of the partial half-lives of isotopes with Z = 117-122. IMELDM is extracted from Ref. [18], 
GLDM+RHF, KPS, and Xu are extracted from [19], Sarriguren is extracted from Ref. [20], Sridhar is extracted from Ref. [73]. 
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fect, as shown in Fig. 2(d, f), compared to £2. Fig. 3 demon- 
strates that the FB increases when Z is large; this indicates 
the competition between the FB and Coulomb repulsion in 
superheavy nuclides. 


The extrapolation stops at the single-neutron (proton) and 
two-neutron (two-proton) drip lines. The UNEDFO data set 
stops at Z — 120. From the existing region to the neutron- 
deficient side, the a decay and SF are predicted to compete. 
On the neutron-rich side, the calculations predict the G~ de- 
cay as the dominant mode, whereas the SF competes for spe- 
cific nuclides. The latest results of most theoretical calcu- 
lations of the partial half-lives [17—20, 22, 23, 70—74] indi- 
cate that the œ decay mode is dominant for new elements at 
N < 184. As shown in Fig. 4, the partial half-lives of iso- 
topes with Z = 117-122 were predicted in this study and 
compared to the corresponding results in Refs. [18-20]. Al- 
though the partial half-lives of 8* decay and EC determined 
in this study were not longer than those indicated in Ref. [20], 
they remain approximately five orders of magnitude greater 
than that of o decay in this region, which does not change the 
dominant decay mode. 
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Fig. 5. (Color online) Odd-even staggering of Tı /2,sF of U, Pu, 


Cm, Cf, Fm, Rf, Sg, Ds, Cn, and Db isotopes and comparison with 
Tija æ 


The T! /2,a values predicted in this study were longer than 
those indicated in Refs. [18-20]; this does not change the 
dominant decay mode of odd-Z isotopes but enhances the 
competition of SF in even-Z isotopes. Furthermore, the pre- 
diction in this study demonstrated a strong odd-even stagger- 
ing of T, 5 sp of even-Z isotopes. In other words, the T; /2 sr 
value of even-even nuclei is several times shorter than its two 
isotopic neighbors, which differs from the weak or unpre- 
dicted odd-even staggering effect obtained by other SF mod- 
els shown in Fig. 4. Notably, all measured T)/2 sr values 
of the even-Z isotopes demonstrate such odd-even stagger- 
ing. Fig. 5 illustrates T1/5.sp and T)/2,, of isotopes with 
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384 


Z > 92. For example, when Z is small in the U, Pu, Cm, 
and Cf isotopes, SF is not competitive with the o decay be- 
cause the Coulomb repulsion is not sufficiently strong. How- 
ever, when Z is large, the Coulomb repulsion increases, and 
this odd—even staggering makes the SF competitive with the a 
decay in the even-even nuclides. Thus, o decay is suggested 
to be a key signal detected for Z — 119 and 121 isotopes, 
whereas the SF should also be considered for even-N iso- 
topes of Z — 120 and 122. Moreover, odd-even staggering 
also exists in odd-Z isotopes, which can only be verified by 
260—263Db because the data are limited. Therefore, the odd— 
even staggering of odd-Z isotopes was not predicted in this 
study. The DNS model predicted the cg value of hundreds 
of FB for the 3n or 2n channels producing ?9??119;74 on the 
243 Am target [75], which can be examined for the new fa- 
cilities of CAFE2 and SHANS2 in Lanzhou [76]. Given the 
odd-even effect of partial half-lives, nuclide candidates for 
new superheavy elements still require an analysis based on 
the cross-section and partial half-life. 


IV. SUMMARY 


In this study, the decay modes of superheavy nuclei were 
investigated using the RF algorithm. The partial half-lives of 
a decay, 8^ decay, 8^ decay, as well as EC, and SF were 
studied and compared. The dominance of o decay in the 
neutron-deficient region was relatively evident. 87 decay is 
predicted to be dominant in the neutron-rich regions. The 
SF contributes to a long-lived circle at the southwest corner 
of Z — 114 and N — 184. A more accurate and precise 
measurement of the nuclear mass and decay energy can im- 
prove the prediction of the decay mode. The odd-even effect 
of the SF was observed in even-Z nuclides. Combined with 
the strong Coulomb repulsion, the SF and o decay became 
competitive in even-even nuclides. Thus, the o decay is sug- 
gested to be a key probe of isotopes with Z — 119 and 121, 
whereas the competition of SF should be considered in even- 
even isotopes with Z — 120 and 122. 

250,252,254 Cy) | 260,261 Es, 261-264) fq, and 2657 7 with half- 
lives predicted to be longer than 10* s were suggested for fu- 
ture measurements. The SF, influenced by the fission barrier 
and Coulomb repulsion, leads to a long-lived region during 
extrapolation. The results of this study indicate that research 
regarding SF, especially beyond ?99FI, which is currently the 
heaviest nuclide with a significant SF branch ratio, is critical 
for performing studies on new facilities, such as CAFE2 and 
SHANS?2 in Lanzhou. 
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